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CN \ ABSTRACT 

The UKIRT Infrared Deep Sky Survey (UKIDSS) is the first of a new generation of 
hemispheric imaging projects to extend the work of the Two Micron All Sky Survey 
(2MASS) by reaching three magnitudes deeper in YJHK imaging, to K = 18.2 (5a, 
Vega) over wide fields. Better complementing existing optical surveys such as the 
Sloan Digital Sky Survey (SDSS), the resulting public imaging catalogues provide new 
photometry of rare object samples too faint to be reached previously. The first data 
release of UKIDSS has already surpassed 2MASS in terms of photons gathered, and 
\Q [ using this new dataset we examine the near-infrared properties of 2837 quasars found 

. in the SDSS and newly catalogued by the UKIDSS in ~189deg 2 . The matched quasars 

include the RA range 22^ to 4' 1 on the Southern Equatorial Stripe (SDSS Stripe 82), 
an area of significant future followup possibilities with deeper surveys and pointed 
observations. The sample covers the redshift and absolute magnitude ranges 0.08 < 
z < 5.03 and —29.5 < Mi < —22.0, and 98 per cent of SDSS quasars have matching 
UKIDSS data. We discuss the photometry, astrometry, and various colour properties 
of the quasars. We also examine the effectiveness of quasar/star separation using the 
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near-infrared passbands. The combination of SDSS ugriz photometry with the YJHK 
near-infrared photometry from UKIDSS over large areas of sky has enormous potential 
for advancing our understanding of the quasar population. 
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1 INTRODUCTION 

In the several years since the SDSS (Sloan Digital Sky Sur- 
vey) and 2MASS (Two Micron All Sky Survey) projects pio- 
neered large uniform digital datasets accessible to the astro- 
nomical community, another generation of surveys has been 
spurred by improvements in imaging technology, software 
pipelines, and the ability to carry out sophisticated astro- 
nomical science on databases. In particular, over the next 
several years, two large consortia will be exploring deeply 
into the near-infrared with surveys from 4-metre class tele- 
scopes. Using UKIDS S (UKIRT Infrared Deep Sky Survey, 
lLawrence et al.ll200fl ) and VIS TA (Visible and I nfrare d Sur- 
vey Telescope for Astronomy, lEmerson et al.l I2004F ) , they 
aim to extend the limiting flux of the 2MASS JHK bands 
by three magnitudes or more, as well as into new parameter 
spaces using the nov el Z and Y optical/near-infrared filters 
|Hewett et al.ll2006l ). 

In part motivated by the search for objects such as the 
coolest brown dwarfs and highest redshift quasars, deeper 



surveys in the near-infrared provide valuable information for 
the practical selection and identification of rare objects from 
parent catalogues of normal stars and other contaminants. 
Both high-redshift quasars and brown dwarfs are marked 
by rising flux towards and beyond 1 fim, but, until recently, 
were rarely discovered in the near-infrared due to the shallow 
limiting magnitude of existing surveys. For normal Galactic 
stars detected in SDSS, for example, 2MASS data has pro- 
vi ded sufficient depth fo r successful characterization, such as 
in lFinlator et al 1 (|2000l ). but brown dwarfs and high-redshift 
quasars are typically much fainter than the 2MASS survey 
limits - only ~17per cent of SDSS-discovered quasars are 
also detected in 2MASS, and these are primarily at low red- 
shifts. 

With the first large data release by the UKIDSS, it is 
now possible to study the near-infrared properties of a large 
sample of SDSS spectroscopically identified quasars over an 
extended redshift range. Further, with an eye towards future 
selection techniques, the addition of deeper near-infrared 
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Figure 1. Positions of SDSS DR3 quasars matched with 
UKIDSS catalog sources are plotted in right ascension (RA) 
and declination (Dec)(J2000). The area plotted is the south- 
ern equatorial stripe (SDSS Stripe 82), a representative sub- 
set of the total area matched in this work. Some obvious 
gaps in RA/Dec, around RA=50° for example, are due to 
currently incomplete UKIDSS imaging, and the over-/under- 
density of matched quasars versus RA (~ 30%) is due to 
the inhomogeneous spectroscopic plate coverage of the input 
SDSS quasar catalogue. (Full-resolution figures will appear in 
published version, or may be downloaded in the meantime at 
http:/ '/www. astro. ex. ac.uk/people/chiu/chiu.figs. tar.gz ) 



photometry has potential for the identification of new classes 
of rare quasars. Just as importantly, deeper near-infrared 
data facilitates the rejection of spurious candidates in such 
surveys. Therefore, combining data from the UKIDSS with 
existing data from SDSS and other surveys may offer large 
gains in observing efficiency. 

In this work, we match spectroscopically identified 
quasars from the SDSS DR3 with their newly released 
UKIDSS DR1 photometry. We discuss the data matching 
procedures in the following sections, examine the joint prop- 
erties of the quasars compared with normal stars, the pho- 
tometric behaviour of the quasars versus redshift, and their 
relation to syn t hetica lly generated model quasars. As in 
ISchneider etafl l|2005h . all absolute magnitude values in this 
paper are based on the standard 'concordance' cosmology of 
n M = 0.3, Q A = 0.7, H = 70 km s" 1 Mpc" 1 . 



2 DATA SOURCES 

2.1 SDSS DR3 Quasar Catalog 

The quasars used in this work are drawn from the SDS S 
DR3 Quasar Catalog, assembled bv lSchneider et all l|2005h . 
The DR3 quasar catalogue was generated from spectroscop- 
ically confirmed quasars discovered and/or reidentified by 
the SDSS over 3732 deg 2 , within the main SDSS northern 
area as well as southern equatorial stripe (SDSS Stripe 82). 
While futher updates to the SDSS quasar catalogue are on- 
going, the DR3 version represents the final product in the 
right ascension (RA) and declination (Dec) regions used in 
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Figure 2. SDSS-UKIDSS matched quasars (bottom panel) are 
shown compared to a sample of raw unidentified point sources 
(top panel) that have also been matched with UKIDSS. Their 
apparent SDSS i and UKIDSS K magnitudes arc plotted, il- 
lustrating the SDSS spectroscopic survey limit (i = 20.2, ver- 
tical dotted line) and UKIDSS photometric limit (K = 18.2, 
horizontal dotted line). (Full-resolution figures will appear in 
published version, or may be downloaded in the meantime at 
http: //www. astro. ex. ac.uk/people/chiu/chiu.figs. tar.gz ) 



this paper. The SDSS automated spectroscopic quasar pro- 
gram, w hose selection criteria and strategy have been dis- 
cussed in iRichards etHl l|2002l ). targets quasars at redshifts 
up to z = 5.4, and down to magnitudes rm < 20.2 (AB), us- 
ing several magnitude and colour cuts applied to the SDSS 
ugriz photometry. 

Following inspection of the individual spectra and com- 
pilation with attributes such as radio, near-infrared, and 
X-ray emission properties, Schneider et al. retained 46,420 
objects, ranging in redshifts from z = 0.08 to 5.41, and 
absolute magnitudes Mi — —22 to —30.2 (including some 
serendipitously discovered quasars down to a spectroscopic 
survey limit of i < 21.0). Typically, the photometric ac- 
curacy of the SDSS is 0.03 mag and systematic astro- 
metric errors are less than 0'.'2. These photometric and 
astrometric properties reflect the high quality of SDSS 
data. Further de tails of the SPSS hardware and s oftware 
can be found in: lYork et al. (2000), Fukug ita et al. 
Gunn et all (| l99 8h,lHogg et all (120011), iLupton et al.1 ( 



Stoughton et alj (I2002TI ISmith et alj i|2002l ), 
l|2003 ). lAbazaiian et al.l (|2004l ). 
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Figure 3. The redshift and absolute Mi magnitude distribu- 
tion of UKIDSS-SDSS matched quasars is shown, and compared 
against the parent SDSS DR3 quasar catalogue. In the central 
panel are plotted the 2837 quasars from the successfully matched 
UKIDSS-SDSS subsample. The slight density edge perceptible 
within the sample results from the magnitude limits of the SDSS 
main quasar spectroscopic survey l|Richards et al]|2002f l. The side 
panels are histograms of the distributions in redshift and absolute 
magnitude, with the parent DR3 quasar catalog shown as a dot- 
ted line, and matched quasars as a solid line. The parent sample 
histograms have been scaled by a factor of 0.4 for comparison. 
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Figure 4. Differences between SDSS and UKIDSS measured po- 
sitions in arcseconds. The dotted line shows results for a bright 
stellar reference population of 34,000 sources in the matched area. 
Solid line shows the sample of 2837 matched quasars. The mod- 
eled profiles have Gaussian FWHM of 0'.'3 for the quasar sample, 
and 0'.'35 for the stellar reference sample. The stellar population 
histograms have been scaled by a factor of 0.14 for comparison. 



2.2 UKIDSS 

UKIDSS represents the first wide and deep ground-based 
survey in the near-infrared regime since 2MASS operations 
concluded in 2001. The UKIDSS0 is composed of several 
sub-surveys targeting different parameter spaces of both 
Galactic and extragalactic observations, planned in a typ- 
ical 'wedding cake' strategy which balances wide-field ver- 



sus deep coverage. The five sub-surveys include deep/narrow 
target fields such as the UDS (Ultra Deep Survey), and 
wide/shallow fields such as the GPS (Galactic Plane Sur- 
vey). Of relevance to this work is the UKIDSS Large Area 
Survey (LAS), which aims to cover an area of 4,000 deg 2 
matching the SDSS over seven years. 

During the ~ 130 nights per y ear in which the WF- 
CAM instrument l|Casali et alj|200lh is mounted on UKIRT 
and devoted to UKIDSS, the project observes the sky un- 
der tolerable seeing conditions (< l'.'2), and photometric 
or mildly non-photometric skies. Each pointing of the tele- 
scope covers an effective area of 0.21 deg 2 , but due to the 
sparse- filled arrangement of the four 2Kx2K infrared detec- 
tors, four pointings are necessary to fill in a full sky patch of 
~ 0.77 deg 2 . Because of the large pixel scale of the detectors 
(O'.'4/pixel) and the large field of view, survey progress is 
quick, covering this filled tile area in the four LAS YJHK 
filters to a depth of K = 18.2 (5a, Vega) in approximately 
20 minutes, with a total of 40 s effective exposure time per 
band. Survey progress under good conditions can thus pro- 
duce appro ximately 2 5 deg 2 of complete LAS i magin g data 
per night. iDve et all (120061 ). lLawrence et al.1 {2006), and 
IWarren et all 12006) provide comprehensive descriptions of 
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the survey hardware, software strategy and data product 
characteristics. 

Following observation, pipeline processing, and cata- 
logue generation, the typical limiting magnitude of the LAS 
is YJHK = [20.16, 19.56, 18.81, 18.19] (Vega, 5a) or [20.79, 
20.49, 20.19, 20.09] (AB, 5a), with photometry having been 
bootstrap-calibrated to the 2MASS system. These limits 
represent an approximately three magnitude depth gain 
compared to the 2MASS JHK imagery. Astrometric accu- 
racy in the LAS is typically < O'.'l, and photometric accuracy 
is ~ 0.04 mag, measured in the J-band. Imaging in the four 
filters is taken in pairs, YJ and HK, for any particular field. 
The imaging catalogues present essentially an instantaneous 
single epoch measurement for each filter-pair but significant 
time intervals can exist between the two filter-pairs. 



3 CATALOGUE MATCHED DATA 
3.1 Matching Procedure 

The UKIDSS DR1 release overlaps a subset of the SDSS 
northern and southern areas and achieves a limiting magni- 
tude nearer to the SDSS for a typical quasar spectral energy 
distribution (SED) than did 2MASS. At the same time, the 
photometric and astrometric performance of the UKIDSS is 
similar to the SDSS. 

The areal extent of the newly released UKIDSS imaging 
is much smaller than the available SDSS DR3 coverage and 
we first trimmed the Schneider et al. (2005) quasar catalogue 
to dimensions similar to the UKIDSS DR1 imaging areas. 
Among the total UKIDSS DR1 imaging area released, a sig- 
nificant portion is located in the southern equatorial stripe, 
which is also a special SDSS region where many repeated 
scans have been taken over the lifetime of the survey. In 
the future, this area may become a target of deeper survey 
programs at various wavelengths. The SDSS and UKIDSS 
have a common southern equatorial stripe area with an RA- 
range of approximately 22 h to 4 h , and a declination-range 
of ±1.265°, centered on 0°. 
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Figure 5. Photometry differences between UKfDSS and 2MASS measurements. Left panel: for a sample of reliable stellar sources 
detected by both UKIDSS and 2MASS with J > 14, the photometry differences between the two surveys are shown versus UKIDSS 
magnitude, in the three 2MASS bands, JHK. Dotted vertical lines in both panels indicate the 2MASS 5cr limiting point-source magnitude 
in each band (JHK = 16.55, 15.85, 15.05, Vega). Right panel: the photometry of 474 (out of 2837 matched) SDSS DR3 quasars found 
in both UKIDSS and 2MASS catalogues is shown, illustrating the sample scatter and that most SDSS DR3 quasars will be found 
below the 2MASS catalogue limits. (Full-resolution figures will appear in published version, or may be downloaded in the meantime at 
| http://www. astro, ex. ac.uk/people/ chiu/ chiu. figs. tar.gz\ ) 



Extracting quasars from the DR3 catalogue without 
any imposed constraints, the quasars were then position- 
ally cross-matched with the UKIDSS DR1 database, with 
a tolerance of l'.'O, matching only the nearest object to the 
queried position. A total of 2837 SDSS quasars were succes- 
fully matched to the UKIDSS DR1 database. 



3.2 Matched sample properties 

3.2.1 General properties 

A portion of the southern equatorial RA/Dec distribution 
of the resulting sample of SDSS-UKIDSS quasars is shown 
in Figure [1] Other sky-areas contributing to the SDSS- 
UKIDSS s ample have more complicated geometries (de- 
scribed in IWarren et aTl 120061 '). but show similar matched 
object densities (within ~ 30 per cent due to SDSS spec- 
troscopic plate coverage over-/under-densities). Within the 
RA-range l h to 2 h 40 m , of 574 input quasars in the south- 
ern equatorial stripe, (corresponding to a sky density of 
9 deg -2 ), 567 (98.7 per cent) have a corresponding UKIDSS 
source with at least one detection in the Y JHK-bands. We 
employ the results of all the matches to the UKIDSS DR1 
catalogue, irrespective of the number of detections among 
the four passbands. Figure [2] illustrates the relationship of 
typical matched quasars to the parent sample of all matched 
point sources in the joint SDSS-UKIDSS dataset, showing 
the spectroscopic and photometric sample limits of the cata- 
logues from which the quasars are drawn. Ongoing projects, 



for example emplo ying the novel AAOm ega multifibre spec- 
trograph at A AT l|Saunders et al.l 12004). aim to extend the 
SDSS-selected quasar sample to magnitudes below these lim- 
its, and improve constraints on the faint and distant quasar 
luminosity function. 

Figure [3] illustrates the redshift and absolute magnitude 
dependence of the full sample of SDSS-UKIDSS quasars, 
which show no strong bias from the parent unmatched sam- 
ple of SDSS quasars. Further, in UKIDSS areas which are 
imaged in all four UKIDSS Y JHK bands, we examined 
the frequency of quasar detection versus band and found 
no strong dependence - when a quasar is detected, it can 
generally be detected in all four bands. Specifically, for ex- 
ample in the above representative area, out of 330 quasars 
found in the smaller 4-band imaged area, 329 are detected 
in the Y-band, 328 in the J-band, 321 in the ff-band and 
325 in the if-band, while 319 have detections in all Y JHK 
bands. 

3.2.2 Morphology 

In the DR3 quasar catalogue, quasars at high redshifts 
are generally un resolved by the SDSS imaging pipeline 
(|Chiu et al.H2005l ). The SDSS-UKIDSS sample shows a sim- 
ilar behaviour in the UKIDSS imaging bands, though with a 
greater fraction of resolved quasars detected at low redshift - 
222 of 278 quasars with redshifts z < 0.5 are classified as ex- 
tended in the UKIDSS mergedClass parameter. Because the 
mergedClass parameter combines weighted measurements 
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Figure 6. UKIDSS-detected but 2MASS-undetected quasars 
can be recovered in composite from 2MASS imaging. 60 SDSS- 
UKIDSS quasars with Jtjkidss > 18 were selected and their 
2MASS atlas images in JHK retrived. Coaddition of the 2MASS- 
images yields the composite source above. Each image is 1' square. 

of the measured source profile from each of the available 
bands, contribution from low redshift quasar host galaxies, 
particularly in the K band, may influence this trend (see 
Maddox & Hewett 2006). However, by z > 1.0, only ~12per 
cent of quasars are resolved in any UKIDSS band, falling to 
7 per cent and 5 per cent by z > 2.0 and z > 3.0 respec- 
tively. In the present sample (without regard to redshift), the 
majority of the quasars are classified as point sources with 
good point spread function (PSF) fits (UKIDSS database 
mergedClass = — 1, N = 1675/2837), many are classified 
as extended (mergedClass = +1, N — 717/2837), and a 
smaller fraction possess poor PSF-fits (mergedClass = —2, 
TV = 435/2837). Ten objects appear as uncertain extended- 
source fits (mergedClass = —3), or are classified as noise. 
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Figure 7. Top panel: throughputs of the SDSS and UKIDSS 
imaging systems. Note the novel UKIDSS Z and Y passbands, 
which have sharper cutoffs compared to SDSS z, and fill the wave- 
length gap between the SDSS and near-infrared JHK bands. For 
reference, the top axis scale indicates the wavelength position 
of quasar Lya emission versus redshift. Bottom panel: horizontal 
lines indicate the resulting 5<r F\ sensitivities (erg/s/cm 2 /A) and 
wavelength coverages of the SDSS and UKIDSS survey bands, in 
addition to the shallower 2MASS in JHK . A composite quasar 
spectrum shifted to redshift 2 = 3 and normalized to magnitude 
iab = 19.1 is shown, demonstrating near-infrared detection in 
UKIDSS but not 2MASS. 



3.2.3 Astrometry 

As cosmologically distant point sources, quasars should 
provide a fixed reference frame on the sky. Therefore, com- 
parisons of the astrometry between different surveys of these 
same objects can yield a measure of the relative astromet- 
ric errors. Using the sample of SDSS-UKIDSS quasars then, 
the positional offsets between the surveys was determined. 
The resulting distribution of relative separations is shown 
in Figure 2] along with a large reference sample of reliably 
measured normal Galactic stars (N ~ 30, 000). While a max- 
imum offset of 0'.'5 can be seen for a few quasars (and up 
to 0'.'6 for the stellar population), the bulk of the sample is 
well modelled by a Gaussian distribution with a full-width 
at half-maximum (FWHM) of 0'.'3. In comparison, the stel- 
lar population displays a somewhat wider distribution of 
astrometric residuals, with a FWHM of ~ 0'.'35. These pa- 
rameters are c onsistent with the SD SS/2MASS astrometric 
comparison m iFinlator etall (120001 ). We also find a small 
systematic absolute offset in both RA and Dec of ~ 0'.'05 
between the SDSS and UKIDSS positions. 

3.2.4 2MASS magnitudes 

In addition to providing SDSS-measured properties, the 
DR3 quasar catalogue was matched with external datasets, 
such as FIRST and ROSAT detections, Galactic extinction 
measurements and 2MASS photometry. Using this infor- 
mation, we examined the differences in quasar photometry 
between the 2MASS and UKIDSS measurements. For the 
SDSS-UKIDSS quasars, we extracted any available 2MASS 
magnitudes in the JHK bands. As can be inferred from the 



resulting sample in Figure [5] many SDSS-UKIDSS quasars 
fall well below the nominal 2MASS 5<r limiting magnitudes 
(JHK = 16.55, 15.85, 15.05, Vega). However, a small frac- 
tion are still measured in the 2MASS catalogues, down to 
approximately J ~ 17.8, H ~ 17.2, and K ~ 16.5. In 
our sample, 474 of 2837 SDSS/UKIDSS matched quasars 
(17 per cent) were sufficiently bright to be detected in the 
2MASS catalogue in any band. Comparing UKIDSS and 
2MASS magnitudes in this range then, we find a scatter of 
a — 0.3 mag between the two quasar measurements (above 
the 5a 2MASS limits), and similar to the envelope of scatter 
in reference point sources at these faint magnitudes, which 
for example ranges from a = 0.1 mag at the bright end 
to a = 0.4 mag at the 5a 2MASS detection limits (Figure 
[5] left panel). While quasar variability may account for a 
small portion of the p hotometric differences, on the order 
of Am ~ 0.15 mag, as IVanden Berk et all (|2004h have dis- 
cussed previously, here the photometric errors should dom- 
inate any such variability contribution. 

Although many quasars are individually not detected in 
the 2MASS imaging, the UKIDSS photometry now allows us 
to demonstrate the latent detection ability of 2MASS. Se- 
lecting 60 SDSS-UKIDSS quasars at J > 18 (which cannot 
be detected in 2MASS), we extracted the 2MASS JHK at- 
las images around each SDSS-UKIDSS quasar position. The 
2MASS atlas images were then stacked, yielding the compos- 
ite source shown in Figured] In areas of sky that will not be 
covered by UKIDSS, or other future surveys, this technique 
may allow 2MASS to provide useful information on popula- 
tions of faint sour ces (a technique recently demonstrated by 
I White et ai]|2006l for example). 
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4 QUASAR COLOURS 

Using the SDSS-UKIDSS quasars, we now examine 
optical and near-infrared colour properties of the sam 
their relation to stellar properties in the same paramc 
space, and consider the utility of quasar-star separation, 
important preliminary note here is that in all calculatk 
as well as resulting figures and tables in this work, we 
the magnitude systems returned automatically and with 
modification by the SDSS and UKIDSS survey databases 
the case of the SDSS, these are asinh magnitudes on the 
system (specifically psfMags), while for the UKIDSS, tf 
are Vega magnitudes calibrated to 2MASS (YJHKaperm 
for point sources). We adopt this approach for the east 
observers querying the joint dataselQ so that AB-to-V 
magnitude and colour transformations (which are sometii 
accidentally neglected) are not necessary for immediate 
terpretation of retrieved data in conjunction with resi 
shown here. Also, slight variations in the choice of convert 
factors (depending on the spectral template used) will thi 
fore not affect the results discussed. For subsequent analj 
however, a single more consistent physical AB magnit 
scheme may be preferable , using AB-Vega tra nsformati 
such as those tabulated in lHewett et alj l|2006l ). 
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4.1 The UKIDSS filter system 

The filter systems of the SDSS and UKIDSS are shown in 
Figure [7] In particular, while the SDSS ugriz and MKO- 
based JHK filters are well-known and commonly used, the 
addition of the novel Z and Y filters of the UKIDSS are 
a r elatively recent de v elopm ent and noteworthy. Proposed 
by IWarren fc Hewettl (|2002T ), and discussed extensively in 
iHewett et al.l ^OOrj ). the UKIDSS Z and Y filters improved 
upon the existing broadband coverage between the optical 
and near-infrared in several respects. 

The Z and Y bands are defined by interference filters 
with sharp wavelength cutoffs and, coupled with the high 
sensitivity of modern infrared detectors, remedy some of 
the deficiencies in current imaging around ~ f/im. First, 
the UKIDSS Z filter is much more sharply defined than 
the SDSS z, because the SDSS z filter itself in fact has a 
relatively flat transmissio n of ~ 0.4 which cont inues out be- 
yond f/im (illustrated bv lFukugita et al"1ll996T ). The overall 
z-band redward cutoff (i.e. total system throughput) then re- 
lies on the declining CCD silicon sensitivity to terminate the 
band, and thus results in a very extended filter profile. The 
CCD-detector sensitivity envelope, beginning with a min- 
imum in the ultraviolet, rising to a maximum around the 
r-band, and terminating with a minimum at 1/xm, largely 
accounts for the highly peaked shape of the SDSS through- 
put. This is in contrast to the relatively flat response of the 
UKIDSS system, owing to the roughly constant sensitivity 
of the HAWAII-2 HgCdTe detectors throughout the near- 
infrared. 

Continuing redward, the F-band (used by the UKIDSS 
LAS) provides a new region of sensitivity in the AA ~ 1500 
A wavelength gap between the z(or Z)-band and the MKO 
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Figure 8. Colours of 2837 SDSS-UKIDSS quasars versus 
redshift, in ugrizY J H K combinations of interest. The pan- 
els generally progress from optical to near-infrared colours. 
The colours are calculated using the natural magnitude sys- 
tems of each survey, without conversion between AB or Vega 
(see 2J, i.e. SDSS magnitudes are AB (asinh) magnitudes 
while UKIDSS magnitudes are Vega-based. The colour-redshift 
tracks of a composite quasar are overlaid (solid line), generated 
from a modified v e rsion of the model quasar SED descibed in 
iMaddox fc Hewettl (2006). (Full -resolution figures will appear in 
published version, or may be downloaded in the meantime at 
http: //www. astro, ex. ac.uk/people/chiu/chiu.figs. tar.gz ) 



J band. Among other motivations, the Y filter was specifi- 
cally designed with the selection of very cool brown dwarfs 
and high-redshift quasars in mind. The increasingly red flux 
of these objects versus temperature (for brown dwarfs) or 
redshift (for quasars) has prove n highly useful in r are object 
searches in the last fe w years (|Chiu et alJl200d iFan et al.l 
120061 : IChiu et alj|2005h , but has reached a barrier as current 
surveys, such as the SDSS, no longer provide sufficient sensi- 
tivity in the 1/im range. The high throughput of the UKIDSS 
system at these wavelengths, combined with sharper band 
profiles, allow the more precise discrimination of potential 
quasar candidate redshifts and brown dwarf types, while 
avoiding unnecessary sky background contamination. 

The resulting F\ sensitivity of the UKIDSS imaging 
is shown in the bottom panel of Figure [7] and compared 
with the SDSS as well as 2MASS sensitivities in the JHK 
bands. A composite quasar spectrum at z = 3 is also 
shown, illustrating the important capability of UKIDSS to 
detect the typical quasar near-infrared continuum, versus 
non-detection in 2MASS. 
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Figure 9. The colour differences between 2MASS-detected 
('bright') and UKIDSS-detected ('faint', 2MASS-undetected) 
quasars are plotted, versus redshift in bins of size Az = 0.2. 
2MASS-detected quasars are found up to z = 3.8, providing the 
upper redshift limit to the comparison sample. Errors are cal- 
culated from the quadrature addition of the standard deviations 
of the colour distribution in each 2MASS and UKIDSS bin. In 
the J — K panel, two histograms (faint lines) show the relative 
distribution of 2MASS and UKIDSS quasars contributing to the 
comparison samples. Typical numbers of quasars in the 2MASS 
sample bins are N = 50 at z = 0, decreasing to N = 5 at z = 3.8, 
while UKIDSS-only detected quasars peak around at Z. = 1.4 
with N = 130. Bins with insufficient numbers for comparison are 
omitted, such as at z = 3.0. 



4.2 Quasar colours vs. redshift 

In Figure we plot the 2837 SDSS-UKIDSS quasars, in 
various colours versus SDSS spectroscopic redshift. In each 
colour, only those quasars with detected flux in both of the 
two contributing bands are plotted, i.e. objects with non- 
detections in the necessary filters (indicated as such as by 
mag = -9999.99 in the SDSS or UKIDSS catalogues) are 
omitted. 

As has been discussed in iRichards et all l|2002l) . 



Chiu et all (120051), iMaddox fc HewettJ (|2006l ), and 
Richards et all (|2006aT ). the behaviour of quasar colours 



is suprisingly well-modeled across extended redshift and 
luminosity ranges by a single composite SED. The passage 
of the strong emission lines in and out of the various filters 
controls the small-scale fluctuation of the quasar colours 
versus redshift, while the extended underlying continuum 
provides the distinguishing characteristic separating quasar 
colours from those of stars. Importantly, for certain ap- 
plications , such as the z > 6 i-dropout selection used by 
iFan et all (|2006l ). absorption by the intervening intergalac- 
tic medium provides a strong spectral discontinuity allowing 
the identification of quasars. This effect can be seen for 
example in the rapidly rising r ~ i colour at redshifts z > 4, 
as more quasar flux in the r band is absorbed by the 
intergalactic medium with increasing redshift. 

The observed quasar sample beyond z ~ 2.5 is af- 
fected by the declining space density of the quasar popu- 
lation and the small volumes probed. As a result, the flux- 



limited sample becomes quite sparse (Figure [8}. In order 
to extend the useful predicted colour-redshift values fur- 
ther than z ~ 2.5, we generated synthetic photometry and 
colours of a composite quasar spectrum to predict the me- 
dian redshift-colour evolution track. The technique has been 
descr i bed previously ('quasar cloning' - IMaddox fc Hewett 
| 2006j ; IRichards et~al]|2006al : iHewett et al.ll2006l ; IChiu et al. 
l2005h . and we follow a similar approach here. The modi- 
fied composite quasar based on the spectrum of Maddox & 
Hewett was employed as the model SED, updated to in- 
clude a base continuum contribution, which dominates at 
A < 10000 A, consisting of a dual power-law form 



-0.5 A < 2600A 
-0.2 A > 2600A, 



(1) 



where a specifies the quasar continuum slope, /„ oc v a . 
A blackbody component, with T= 1850 K, dominates at 
A > 10000 A, resulting in an overall spectral shape simi- 
lar to the composite spectrum of Glikman et al. (2006). For 
each redshift bin in steps of Az = 0.1, the model SED was 
redshifted and the average effect of intervening Lyman-a ab- 
sorption added. The resulting quasar spectra were used to 
generate magnitudes and colours, again with SDSS mag- 
nitudes on the AB system, and UKIDSS magnitudes on 
the Vega system, as would naturally be found in the re- 
spective survey databases. These are overlaid in Figure [5] 
Similar colour-redshift loci, based dire ctly on the various 
composite quasar spect r a available (e.g. Francis et al.|[l99ll ; 
iBrotherton et all |200ll ; IVanden Berk et al.1 1200 ll ). are fre- 
quently employed for a variety of purposes. It should be 
stressed that the significant contribution of host galaxy light 
in the low-luminosity quasars that typically dominate the 
composite spectra at wavelengths > 5000 A leads to predic- 
tions that are far too red in the optical-near-infrared colours, 
compared to the observations, at redshifts z > 2. 

In addition, given the large magnitude range covered by 
the SDSS-UKIDSS quasars, we investigated whether the ob- 
served distribution of quasar colours showed dependence on 
luminosity. This was accomplished by selecting two subsets 
of quasars, those bright enough to be detected in 2MASS, 
and the remaining ones detected only in UKIDSS. The 
colour differences between the two samples were calculated 
in redshift bins of size Az — 0.2 and are shown in Figure [9] 
No significant difference is seen between the two samples. 



4.3 Quasar colours versus normal stars 

To compare the quasar colours with those of normal stars, 
we extracted a stellar sample from the UKIDSS/SDSS 
database. This was accomplished by querying the UKIDSS 
reliableLasPointSource table, which returned relatively 
bright objects displaying good photometric errors (a < 0.1 
mag), unresolved shape/flux profiles, and in uncrowded ar- 
eas. Then, for a few selected two-colour spaces of interest, 
the stellar and quasar samples are plotted in Figure [101 

Photometric searches for quasars generally compete 
against heavy contamination by stellar objects with similar 
colours, due to the intrinsic rarity of quasars. The relative 
distributions of quasars versus stars in a chosen colour com- 
bination will generally control the search efficiency, i.e. how 
many contaminants must be observed to return one object 
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Figure 10. SDSS-UKIDSS quasars (small dots) in colour-colour diagrams (ugrizY J H K) with the density of normal stars shown as 
contours. Quasar-star separation is highly effective in colours involving the bluest bands, such as u — g, but then decreases in utility in 
the optical. With near-infrared and particularly X-band imaging, quasar-star separation again can provide a powerful selection method. 
Note: colours here are calculated from natural magnitude systems of each survey, without conversion between AB or Vega — i.e. SDSS 
magnitudes are AB (asinh) magnitudes while UKIDSS magnitudes are Vega based (see Q. Quasar symbol colours are encoded as follows 
— black: z < 0.5, blue: 0.5 < z < 1.0, green: 1.0 < z < 2.0, yellow: 2.0 < z < 3.0, red: 3.0 < z < 6.0) (Full-resolution figures will appear 
in published version, or may be downloaded in the meantime at\http://www. astro. ex. ac.uk/people/chiu/chiu.figs.tar. gz\) 



from the desired sample. The progression of the two-colour 
plots in Figure [TU] illustrates the varying success of distin- 
guishing quasars from normal stars using different colours. 
For colours involving the ultraviolet, such as in (u — g), 
(g — r), the strong short- wavelength quasar continuum flux 
is well known to provide a separation from stars via the 
ultraviolet-excess (UVX) technique. Relatively clear separa- 
tion of the two populations is observed here and provides a 
highly successful selection method for low redshift quasars. 

In the optical bands, adjacent filter combinations pro- 
vide little discrimination between the quasar and stellar pop- 
ulations. However, with the addition of near-infrared data, 
and particularly K-band imaging, the extended quasar con- 
tinua again results in an effective differentiation between 
quasars and stars. This att ribute, exploited by the 'KX' 
method dWarren et alj|2000h . can provide a useful selection 
at high redshift, z > 2.5, where the UVX technique is un- 
successful due to the effect of absorption by the IGM in the 
bluest passbands. 

Looking farther into the infrared, an additional area 
of significant interest in the future may be the combina- 
tion of deep near-infrared surveys such as UKIDSS with 



mid-infrared legacy imaging from the Spitzer Spa c e Tele - 
scope. As has been demons trated by Stern et al] ll2005h , 
lHatziminaoglou et al.l (|2005h . and iRichards et al] (<2006bh . 
quasar-star separation is highly successful in the mid- 
infrared bands, and also benefits from high sensitivity due to 
the supressed background (15/xJy in 90s at 3.6/im). In areas 
with good optical, near-infrared, and proposed mid-infrared 
imaging coverage from Spitzer, such as the southern equato- 
rial stripe, the combination of the resulting catalogues has 
the potential to reach new distant and faint quasar popula- 
tions while providing nearly continuous bandpass coverage 
as emission lines and continuum flux shift with redshift. 



5 QUASAR COLOUR SELECTION 
STRATEGIES 

The utility of quasars in probing the distant universe has 
been demonstrated to great effect in the past few years, 
with these high redshift objects illuminating not only their 
own properties, but also the physical state of the IGM and 
the distant universe. However, the discovery of quasars in 
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Figure 11. A detailed view of the (J—H), (H—K) distribution of 
stellar sources (contour lines) and SDSS-UKIDSS quasars at 2.2 < 
z < 3.3 (red + symbols), illustrating the useful colour separation 
of populations with i\-band information. A model quasar track 
(discussed in text, blue solid line) is shown from z = 0.0 to z = 5.5 
with Az = 0.1 interval ticks and numeric labels. The quasar track 
overlap with stellar sources is severe at z ~ 1.4, but reaches large 
separation from z = 2.0 to z = 3.5. 



certain redshift intervals is, in practice, a laborious enter- 
prise. In the absence of information from the wavelength 
extremes, identifying quasars by special properties such as 
radio or X-ray flux, discovery of quasars at redshifts near 
z ~ 2.5, z ~ 5.6, and z > 6 remains a significant challenge. 
Observers fight the combination of the intrinsically declining 
numbers of the target population and the increasing num- 
bers of contaminating Galactic stars. In previous work we 
have outlined some of the selection techniques and num- 
bers of quasars to be expected in the SDSS and UKIDSS, 
and have found that the UKIDSS-SDS S joint dataset is ca- 
pable of selecting quasars at z > 6 5 dHewett et al.ll2006l : 
IChiu et al.ll2005l ; IWarren fc HewettH2002t ). However, at the 
magnitude limits of the UKIDSS (m A B ~ 20 - 21), the 
N ~ 90 expected quasars at z > 6.5 on the whole sky 
clearly places this parameter space in the category of rare, 
high-value, but resource-intensive targets. 

Here, we discuss the more tractable question of improv- 
ing the effectiveness of quasar selection at redshifts z ~ 2.5 
and z ~ 5.6 where colour-selec ted samples in the o ptical 
encounter particular difficulties (|Schneider et al.ll2005l ) . The 
relatively deep UKIDSS, and future VISTA, limiting mag- 
nitudes in the near-infrared, and the high rate of success in 
matching the present sample of SDSS quasars to UKIDSS, 
suggests that quasars at z ~ 2.5 and z ~ 5.6 may be recov- 
ered by employing if -band information. Figure[TT]shows, for 
example, a more detailed view of the quasar-star separation 
in colour space when Tf-band information is added. While 
the population overlap (and thus contamination) is signifi- 
cant for predicted quasar colours around z ~ 1.4, the colour 
separation becomes large from z = 2.0 to z = 3.5. Because 
the quasars shown are gathered from optical selection tech- 
niques, this implies that independent (or joint) use of the 
near-infrared photometry on unidentified stellar sources in 
such a region should yield additional quasars. 

Thus, an appropriate choice of optical and near-infrared 
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Figure 12. An example of a successful low-redshift quasar se- 
lection using colour separation in (u — z),(Y — K) space is il- 
lustrated. Normal stars are plotted as small dots, while SDSS- 
UKIDSS quasars are shown as X symbols, with redshift values 
encoded by symbol colour as in Figure HOI A quasar track, from 
z = 0.0 (at top) to z > 3.5 (at far right), with numeric labels 
at intervals of Az = 0.4, derived f rom a modified version of the 
quasar SED of iMaddox &" Hcwct is also shown. 

colours can select quasars in the redshift ranges of interest. 
Figure[l2]illustrates one example of a successful quasar iden- 
tification/separation scheme in the (u — z),(Y — K) colour 
space which more cleanly separates normal Galactic stars 
from quasars by sampling the underlying char acteristic spec- 
tral b reak / inflexion featur e s, as d escribed bv lWarren et al.l 
|2000l ) and lRichards et all (|2006bl l. 

This colour selection boundary illustrated is given by 
the following criteria: 

(u-z)< 4.5 (2) 
(Y-K)> 0.6 (3) 
(Y - K) > 0.35* (u-z) + 0.425, (4) 

which results in a completeness of 97 per cent (measured 
by the fraction of quasars selected by the above cut) for those 
quasars with magnitudes in these bands. The overall success 
of the selection will vary as a function of the magnitudes of 
the target population compared to the limiting magnitudes 
in each passband. For instance, the selection of z > 3 quasars 
would be most sensitive (and efficient) using colour cuts in 
the SDSS g band and beyond, such as (r — Y), (Y — K). 

As illustrated above, the availability of photometry cov- 
ering the rest-frame optical and near-infrared wavelength re- 
gions offers the prospect of obtaining a more complete census 
of the quasar population, both in terms of the redshift range 
accessible and the variety of quasar SEDs included. A conse- 
quence of the extended restframe wavelength range accessed 
through obervations spanning the u through K bands for 
quasars with redshifts < z < 7, is the importance of incor- 
porating realistic model quasar SEDs in order to interpret 
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UKIDSS DR1 source catalogue. The resulting sample of 
2837 quasars is representative of the full range of magni- 
tudes and redshifts in the SDSS spectroscopic sample. The 
astrometric properties of the sample are similar to previous 
experience with 2MASS-SDSS matched sources, while the 
UKIDSS photometry extends the matched imaging depth by 
three magnitudes in the JHK bands compared to 2MASS. 
The median colours of the matched quasars agree well with 
predictions generated using a model quasar SED and pro- 
vide a guide for higher redshift selection techniques. Opti- 
cal and near-infrared colours, in appropriate combinations, 
provide effective separation of stars and known quasars over 
essentially the full redshift range < z < 7, overcoming the 
limitations of employing only optical colours in the redshift 
ranges z ~ 2.5 and z ~ 5.6. 



Figure 13. As quasars of increasing redshift are selected us- 
ing the KX-technique, the K filter probes increasingly blue rest- 
frame wavelengths of the quasar SEDs. As a result, different 
continuum slope regions (et pt/iiir) are sampled and it is im- 
portant that appropriate i^-corrections are incorporated in any 
analysis of the properties of the target po pulation. A rest-frame 
qua sar composite spect r um d erived from iGlikman et al ] l|2006l l 
and IVanden Berk et al.l ll200ll ') is shown to illustrate the impor- 
tance of incorporating information about the SEDs of the target 
population(s). 

the properties of the quasar population reliably. The sim- 
ple power-law approximation to quasar SEDs, modulated 
only by the presence of the strongest emission lines, that 
is relatively effective when interpreting observations prob- 
ing only the restframe wavelength region AA1200 - 5000 A 
of high-luminosity quasars over redshifts 0.5 < z < 2.5 is en- 
tirely indequate when analysing quasar samples of the type 
presented here. To illustrate the point graphically, Figure 
1131 shows a composite quas ar spectrum constructed from 
the ultraviolet composite of IVanden Berk et al.l (|200lh and 
the optical an d near- infrared portion of the composite of 
IGlikman et"aH (|2006l ). joined at 3100 A. The very large 
colour changes in the optical-near-infrared colours as a func- 
tion of redshift evident in Figure [S] are a direct manifes- 
tation of the form of the (very non-power-law-like) quasar 
SED. Further care is required when considering the effect 
of host galaxy contamination on the observed properties of 
the 'quasar s', particularly at redshif ts z < 1 and at lower lu- 
minosities l|Maddox fc Hewettll2006h . The availability of the 
deep near-infrared photometry would also allow a quasar 
sample to be selected over a broader redshift range using 
observed passbands that closely approximate a fixed rest- 
frame wavelength range for the quasars at each redshift, 
thereby minimising the k-correction factors. However, the 
more ambitious goal is now to utilise the availability of the 
full set of ugrizYJHK bands to define samples that will 
provide powerful new constraints on the properties of the 
quasar population as a function of redshift, luminosity and 
SED- type. 

6 CONCLUSION 

We have matched 98 per cent of SDSS DR3 quasars over 
~ 189 deg 2 with their near-infrared counterparts in the 
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